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Abstract 

We propose a new entanglement switch of qubits consisting of electric dipoles, oriented along 
or against an external electric field and coupled by the electric dipole-dipole interaction. The 
pairwise entanglement can be tuned and controlled by the ratio of the Rabi frequency and the 
dipole-dipole coupling strength. Tuning the entanglement can be achieved for one, two and three- 
dimensional arrangements of the qubits. The feasibility of building such an entanglement switch 
is also discussed. 
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Entanglement is a quantum mechanical property that describes a correlation between 
quantum mechanical systems. It has no classical analog and has been lying in the heart 
of the foundation of quantum mechanics. The desire to understand, tune and manipulate 
quantum entanglement is of fundamental importance in the field of quantum information and 



computation 



]. Recently, we studied a set of localized spins coupled through exchange 



interaction and subject to an external magnetic filed 



We demonstrated for such a 



class of one-dimensional magnetic systems, that entanglement can be controlled and tuned 
by varying the anisotropy parameter in the Hamiltonian and by introducing impurities into 
the systems {4]. In this letter, we propose a new entanglement switch of qubits consisting of 
the electric dipole moment of diatomic polar molecules, oriented along or against an external 
electric field and coupled by the electric dipole-dipole interaction. 

Recent progress in methods for producing, trapping and controlling cold polar molecules 
make them an excellent candidate for quantum computation 9Ml2| . Trapped polar molecules 
was proposed as a novel physical realization of a quantum computer by a number of authors 
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-|23[. In this proposal, the qubits are the molecular electric dipole moments which 
can only orient along(|0)) or against (|1)) the external electric fields. Each qubit is one polar 
molecule with equal spacing along the axis. For such a system, tuning and controlling the 
entanglement between the dipoles is of great importance. 

The Hamiltonian of N-trapped dipoles in an external electric field reads 

N N 

n = hY,u l s* + hY,^s+s-, (1) 



i=i 



where S"s are the dipole operators and related to Pauli matrices, and Sj represent dipole 
excitation and de-excitation, respectively. Ui is the transition frequency of the dipole on site 
i, which is a function of dipole moment and external electric fields at site i. 



d-E 



(2) 



where d is the electric dipole moment which we assume the same for each site. fiy is the 
dipole-dipole coupling constant between sites i and j, which is determined by the strength 
of dipole moment, the direction of the external electric field and the lattice constant of the 
dipole array. 

\d\\l 



(3) 



where 6 is the angle between and external electric field. 

We will expand the Hamiltonian in pure standard basis, {|00), 1 01) , |10), |11)}, and all 
the eigenstates will be obtained by diagonalizing the Hamiltonian matrix. For example, for 
iV = 2, we obtained the following four eigenvectors: = |00), |^2 = ( 1 01) — 1 10) ) , |^ 



'3 



4g(|01) + 1 10) ) , \l/4 = 1 11) with corresponding eigenvalues: E\ — 0, E 2 = uj — Q, E 3 = 
uj + Q, E4 = 2u. For the ratio oj/Q < 1, ^2 is fully entangled ground state. However, for 
u/fl > 1 the non entangled state ^1 is the ground state. Thus, we have a curve crossing at 
u/Vt = 1. 

The concept of entanglement of formation is related to the amount of entanglement 
needed to prepare the state p, where p is the density matrix. It was shown by Wootters 



that concurrence is a good measure of entanglement 2QJ]. The concurrence C is given by 



C(p) = max{0, Ai - A 2 - A 3 - A 4 } (4) 

For a general state of two qubits, Aj's are the eigenvalues, in decreasing order, of the Her- 
mitian matrix R = A y ' ^fpp^fp where p is the density matrix and p is the spin-flipped state 
defined as 

p = {a y ®(Ty)p*{<7 y ®(Ty) (5) 

where the p* is the complex conjugate of p and is taken in the standard basis, which for a 
pair of 2 level particles is {|00), |01), |10), |11)}. 

In order to calculate thermal entanglement, we need the temperature dependent density 
matrix and the density matrix for a system in equilibrium at a temperature T reads: p = 
e -l3H/z p — _i_ anc [ 2 is the partition function, Z = Tr^ 1311 ). In this case, the 

partition function is 

Z(T) = J>e-^ (6) 

i 

where Aj is the ith eigenvalue and gi is the degeneracy. And the corresponding density 
matrix can be written: 

1 N 

i 

where |$j) is the ith eigenfunction. For pairwise thermal entanglement, we can get reduced 
density matrix as a function of temperature in the same way, which leads to temperature 
dependent entanglement. 
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In Figure 1 we show the tuning of the pairwise entanglement, measured by concurrence, 
of one-dimensional arrangements of the dipoles as one varies the ratio u/Q at different tem- 
peratures for N = 9 dipoles. Here we took all transition frequencies to be the same, Ui = u> 
and the nearest neighbor dipole-dipole interaction to be the same, = = fi- All 

other dipole-dipole coupling constants f^^ii can be expressed in terms of Q. Thus, we 
have two parameters to vary, the ratio oj/Q and temperature kT. At kT ~ one has a 
constant entanglement over a long ratio u/Q and sharp transitions or jumps to lower values 
at other values of uj/Vt. It is worth mentioning that for u/Q < 0.634, entanglement is only 
between dipole 1 and 2. For 0.634 < oj/Q < 1.14 dipole one becomes entangled also with 
dipole 3 and with other dipoles until we reach oj/Q ~ 1.74, above this value the concurrence 
is zero between all sites. As one increase the temperature, the curve become smoother as 
mixing occur with higher states. Calculations for iV = 3,4, ...,8 gave similar behavior of 
tuning and controlling entanglement as for the case N = 9. To show how the populations 
changes at the transition point, we present in Figure 2 the coefficients of the wave function 
for N = A before and after the transition point u/Q = 0.64. When u/Q < 0.64, the ground 
state wave function is = 0.19| ||tt> ~ 0.51| ||> + 0.45| ||tl) + 0.45| t||t> ~ 0.51 1 t||| 
) + 0. 19 1 till), however, when oj/Q > 0.64, the ground state wave function becomes 
= 0.36| |||t) - 0.61 1 ||tl) + 0.61| |tl|) - 0.36| till)- Figures 3 and 4 show a similar 
phenomena for two and three-dimensional arrangements of the dipoles. The pairwise en- 
tanglement decreases as one increases the dimensionality of the system and the temperature. 



There have been rapid advances made recently in coo 
atomic (Rydberg) 24M26I] and molecular dipoles [id. 
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ing, trapping and manipulating 



29] . For example, a wide variety 



of ground state polar molecules with large electric dipole moments (several Debyes, where 
1 Debye ~ 3.3 x 10~ 30 Cm) have been cooled to ultracold (~ mK or below) regime, some 
even near quantum degeneracy [29]. These developments provide exciting opportunities to 
experimentally realize the entanglement switch described above, which requires tuning uj/Vt 
around 1, and cooling the dipoles to temperatures (T) corresponding to a fraction of the 
dipole-dipole interaction energy (Q). For a dipole moment (p) of few Debyes, and a typical 
experimental electric field (~ 10 5 V/m, required for the dipole moment to actually manifest 
3lj]), dipole-dipole separation (d) on the order of 10 nm is required for u/Q ~ 1. 



itself 
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Such d is much shorter than what can be achieved in typical optical lattices (as envisaged 
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12], which corresponds to the regime of u/Q ^> 1), but can be realized with arrays of 



32 



33] (where 



nanoscale plasmon-enhanced electric/electro-optical traps recently proposed 
extremely tight, few-nanometer confinement and trap frequencies exceeding 100MHz are 

2 

shown to be possible). At such short d, the dipole-dipole interaction f2 = ^ d3 is on the 
order of 0.1 K. Cooling to a small fraction of such a temperature scale is easily within the 
current experimental technology. Furthermore, u/Q can be tuned experimentally by the 
electrical field, and/or by d (eg. by varying the microtrap configuration). 

In summary, we presented a new way to construct an entanglement switch in an optical 
lattice and discussed the experimental feasibility of building such switch. The realization 
of such a scheme will have a profound impact on the implementation of quantum gates 
in quantum computing with trapped polar molecules. The similar idea could in principle 
also be applied to other systems with electric dipoles, such as quantum dots, excitons in 
nanostructurs. 
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FIG. 1: Pairwise concurrence (Eq. 4) of one- dimensional arrangements of the dipoles as one varies 
the ratio u/Q at kT = 10~ 4 , and 10 _1 hQ respectively for N = 9 dipoles. 
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FIG. 2: Pairwise concurrence (Eq. 4) an d coefficients of the wave function of one- dimensional 
arrangements of the dipoles as one varies the ratio u/Q at kT = 10 -4 fiQ, for N = 4 dipoles. 
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FIG. 3: Pairwise concurrence (Eq. 4) for different two-dimensional arrangements of the dipoles as 
one varies the ratio u/Q at kT = 10~ 4 h£l. The external electric field is perpendicular to the 2-d 
plane. 
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FIG. 4: Pairwise concurrence (Eq. 4) of three-dimensional simple cubic arrangements of the dipoles 
as one varies the ratio lo/Q at kT = 10 -4 HQ. Here Q is the dipole-dipole coupling constant between 
dipole 1 and dipole 2, VL = Q.\2- 
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